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ABSTRACT 

A gravitational lens system can be perturbed by "rogue systems" in angu- 
lar proximities but at different distances. A point mass perturbed by another 
point mass can be considered as a large separation approximation of the double 
scattering two point mass (DSTP) lens. The resulting effective lens depends on 
whether the perturber is closer to or farther from the observer than the main 
lens system. The caustic is smaller than that of the large separation binary lens 
when the perturber is the first scatterer; the caustic is similar in size with the 
large separation binary lens when the perturber is the last scatterer. Modelling 
of a gravitational lensing by a galaxy requires extra terms other than constanst 
shear for the perturbers at different redshifts. Double scattering two distributed 
mass (DSTD) lens is considered. The perturbing galaxy behaves as a monopole 
- or a point mass - because the dipole moment of the elliptic mass distribution 
is zero. 
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1. Introduction 



The Galactic bulge is being surveyed for gravitational microlensing in search of mi- 
crolensing planets. The lensing systems are the standard single scattering n-point mass 
lenses, where the bound system may consist of a single star, a multiple star, or a planetary 
system with one or two host stars. The Galactic bulge microlensing probability is ~ 1CT 6 
and the probability for an unbound system to be aligned with the main lensing system can 
be ignored because it is ~ 10~ 12 . The number of stars being surveyed is less than 10 9 . 
The lensing probability is proportional to the Einstein ring radius square, and the Einstein 
ring radius of the lensing toward the Galactic Bulge is characteristically ~ 1 mas. Thus, 
if we consider the probability of a "rogue" system to be within 1 as from the main lensing 
system, the probability is ~ 1. In fact, ground-based microlensing events are known to be 
"plagued" with blending of light, and some of them other than the main lens itself can also 
be gravitationally relevant for the photon path. The perturbers can be dark as well. 

The probability for the "rogue" system to be at the same distance as the main lens 
system is small where the "same distance" should be understood in the context that the two 
lens systems are within the coherence scale in which the lensing can be considered a single 
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scattering lensing. Within the coherence scale, the probability for the photon path to weave 



through the two lenses can be ignored (Rhie and Bennett 2010) (RB10 from here on). Thus 



it is most reasonable to assume in general that the two gravitationally unbound microlens 
systems are at different distances, and they would be best considered as a double scattering 
lensing system. Here it is assumed that the two lens elements are widely separated in the 
sky based on the argument in the previous paragraph and calculate the effects of the "rogue" 
system on the main lens. The double scattering lensing is a time-sequential process, and it 
matters whether the "rogue" system is farther away than the main lens from the observer or 
closer. The two cases are schematically shown in figure [TJ By wide separation it is implied 
that the separation I is much larger than the Einstein ring radius of the main lens. 

We assume that the "rogue" perturbing system is a single point mass. The main lens 
of interest will be a single star, a multiple star, a planet system with one or two host stars, 
or a wide binary stars one of which hosts planets. Here we consider the most common and 
simplest case of a signle star and study the wide separation approximation of the double 
scattering two point mass (DSTP) lens. Then the most important effect of the perturbation 
is to break the degeneracy of the point caustic of the single lens to an extended caustic curve, 
and the size of the caustic will be the indicator of the influence of the perturber. It will be 
shown that the caustic size depends on the distances of the lenses and whether the perturber 
is in the front or in the back. When the perturber is the first scatterer, the caustic is smaller 
than that of the binary lens, and it is similar in size when the main lens is perturbed by a 
"rogue" system in front. A binary lens forms when the two lenses have the same distance - 
or within the coherence length. The binary lens at large separation is made of a point mass 
and a constant shear (and plus the source shift), and it is briefly discussed in the appendix. 

A multiple-point mass lens perturbing a single point mass main lens is approximated 
by the same form of the approximate DSTP lens equation (with effective coefficients) and 
can be concluded to behave in the similar manner to the single point mass perturber. 

In lensing by a galaxy, modelling is done customarily by assuming an elliptic mass 
(sometimes replaced by an elliptic potential) and a constant shear. The galaxy lensing 
is of order 1 arcsecond, and there are often other galaxies in the angular vicinity of the 
main lens system. The perturbers can be group memebers of the main lens or galaxies 
at different distances. If there is a perturbing galaxy at the same distance, its monopole 
will add an external shear as is the case with the binary lens. However, if the perturbing 
masses are at different distances, the perturbation should include deflection terms other 
than the constant external shear. We consider the wide separation approximation of the 
double scattering two distributed mass (DSTD) lens. The other terms depend on the double 
scattering parameter and vanish when the perturbers are at the same distance as the main 
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lens because the double scattering parameter vanishes. The perturbation by an elliptic mass 
galaxy is approximated by the perturbation by a point mass (monopole) because the dipole 
moment of the elliptic mass distribution is zero. The main lens galaxy, assumed to have 
an elliptic mass distribution, has a finite size caustic, and the effect of the perturber is to 
change its shape, size, and position. Even in the simpler case of the main galaxy as the 
monopole-quadrupole lens requires numerical calculations. We leave the perturbations of 
the finite size caustics for future work. 



It should be necessary to point out that Keeton (2003) uses Taylor expansion and 
concludes that the effects of a perturbing mass on the galaxy lens is to add constant conver- 
gence and constant shear irrelevantly of whether the perturber is the first scatterer or the 
last. There may be a problem in the expansion cutoff. In the region of interest around the 
critical curve of the galaxy lens, the first term in the Taylor expansion is small because the 
Jacobian determinant is zero or small and is likely to be smaller than the second order term. 
The second order term is well known for the square root behavior of the lensing near the 
critical curve or caustic crossing. It is not clear whether the Taylor expansion can be used at 
all. We use power expansion around the critical curve of the main lens which is the region 
of interest. 

The DSTD lens equation is an obvious extension of the DSTP lens equation in which 
the delta function integral for the 2-d gravitational field of a point mass is generalized to the 
density function integral for the 2-d gravitational field of the distributed mass. 

The DSTP lens equation is known since 1986 (Blandford and Narayan) and have been 
studied ( Kochanek and Apostolakis|1986 Erdl and Schneider||1993 Werner et al||2008 ) . Here 



the derivation of the DSTP lens equation studied in RB10 is briefed for clarity and conve- 
nience. Instead of using the formula for the time delay and Fermat principle, the well-known 
derivation of the single lens equation from an exact solution of the general relativity, the 
Schwarzschild metric, with the assumption of the linear gravity and small angle approxima- 
tion is used. The Schwarzschild metric is asymptotically flat and the scattering planes can 
be joined easily in the asymptotic regions. The DSTP lens equation is obtained by joining 
two scattering planes with the freedom to rotate. 



2. The DSTP Lens Equation 

The double scattering two point lens equation can be obtained from a diagram shown 
in figure [2] where the linear gravity and small angle approximation are assumed. Since the 
true photon path is three dimensional because of the rotation of the scattering planes with 
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respect to each other, a three-dimensional diagram is needed. But it has been shown in RB10 
that in the linear approximation in small angles, the radial component (in the direction of 
the line of sight) of the impact vector that is generated due to the relative rotation of the 
scattering planes can be ignored because it is of the second order. It is sufficient to express 
the triangular relations of the angles in vectors to account for the relative rotaion between 
the two scattering planes. 

From figure [2] two sets of relations are obtained. 

h = D n {a - 71) + (Ai - D l2 )6<p 2 ; A(«i - P) = -<^i(A - Ai) (1) 

h = D l2 {a - 72) ; D s (a - 5i) = -6cp 2 {D s - D l2 ) (2) 
where the bending (scattering) angles are given by the point mass bending angles. 

? 4GM 1 (-fc 1 ) T 4GM 2 (-b 2 ) 

dipi = ; 6cp 2 = 3 

h b 2 

Mi and M 2 are the masses of the first and second point mass scatterers at the distances Ai 
and A25 and A is the distance to the source; b\ and b 2 are the impact vectors, and b\ = \b\\ 
and b 2 = \b 2 \. The lens equation is obtained from the second equations of eqs. ([IJ and §2^, 

A(« -0) = -hi{D s - Ai) - h*{D, - D l2 ) , (4) 

which is completed by using eq.(|3| and the first equations of eqs. ([I]) and ([2]). 

It is convenient (or our custom) to define a lens plane and use linear variables instead 
of the angular variables. Note that the intermediate image position angle 6l\ was defined by 
projecting the intermediate photon ray back to the sky at the distance of the source. So define 
the lens plane, where the lens equation variables are defined, as the plane at the distance 
of the source and normal to a chosen radial direction. The lens equation is independent of 
the choice of the radial direction because of the linear approximation in small angle. Since 
the lens plane is placed at the distance of the source, the linear variables are A times the 
angular variables. 

Now employ the complex coordinates as usual and let u, z, and Xj : j = 1, 2 denote the 
(2-dimensional) positions (on the lens plane at the distance of the source) of a source, an 
image, and lenses 1 and 2. Then the lens equation in eq.Q can be written in terms of the 
linear variables. 



' E I r E2 , - , 

uj = z-- (5) 

Z l ~ T E2\ Z 2 Z 2 

where = 1,2. Let r Ek j be the (single lens) Einstein ring radius of the lensing 

by object k of object j, and let object refer to the source. r E1 = r E10 , r E2 = r E20 , and r E21 
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are as follows. 



rl 2 = 4GM 2 D 2 § = R% 2 §; D . 2 ^M£±z£A (7 ) 



Dl d2 Dl n _ D l2 {D ll -D, 



r E21 = 4GM 1 D 3 ^- = i?| 21 ^; A3 = -" v ~" ~' 2J 
A2 ^22 Mi 

Dj : j = 1, 2, 3 are the reduced distances, and REkj is the "intrinsic" Einstein ring radius of 

the lensing of object j by object k. The reason why we refer to REkj as the intrinsic Einstein 

ring radius is that the photon rays of the Einstein ring image of the lensing (of object j by 

object k) actually pass through the ring around the lens (k) of radius REkj (accurate within 

the small angle approximation). 

Redefine distances Dj% = Dij and Dj 2 = D s — Dij and define effective masses 

'JJU.-U;'^ 1 : J = l,2. (9) 

Define the Einstein ring radius te of the total effecitve mass, 

r% = r 2 E1 + r| 2 = AGD s (effM 1 + effM 2 ) , (10) 

and the lens equation can be normalized so that the unit distance is te- By substituing u, z 
and Xj in eq.(|5]) by r E uj, r E z and r E Xj respectively, the normalized lens equation is obtained. 

uj = z a~yi ( n ) 

Zi Z 2 

Z 2 

where the effective fractional masses are 

r|i = effMi = Mi = _J_ , . 

1 r| effM 1 + effM 2 M 1 + M 2 /d 1 + e 1 ) 

^r% 1= effM 2 M 2 _ J_ 

62 r| effM 1 + effM 2 M x d + M 2 1 + e 1 j 

and the double scattering parameter is 

_r| 21 M 2 (l-d) 1-d 



The distance parameter <i is 



r 2 , j^d + M 2 1 + e 



(14) 



D 12 D 21 

d=—— < 1 (15) 

D\\L) 22 
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where the equality in the second relation holds when the two lenses are at the distance. The 
effective mass ratio e is 

_ ex _ r% x _ effM x _ Mx 
€ ~e 2 r| 2 e//M 2 V 2 libj 

The effective mass ratio is smaller than the mass ratio. The weight is shifted to the last 
scatterer in double scattering lensing. Note that we have chosen the last scatterer for the 
reference mass. 

It should be worth pointing out that the double scattering parameter a is essentially 
the (square of the) Einstein radius that is easily measurable in an (almost) axisymmetric 



system as in SDSSJ0946+1006 (Gavazzi et al 2008). In an axisymmetric DSTP lens three 



ring images are formed, even though the innermost ring is "unstable" to break into a half- 
circle, and \fa measures the middle ring radius in units of the Einstein ring radius of the 
total effective mass te- The DSTP lens system can be considered to have two characteristic 
parameters Te and \/a. 



Here the focus is in the main lens and the interest is on what happens to the Einsteing 
ring of the main lens under the perturbation of a perturbing mass. So it is useful to renor- 
malize the lens equation by the Einstein ring radius of the main lens. There are two cases: 
1) object 1 is the perturbing mass; 2) object 2 is the perturbing mass. 



Case 1): Renormalize the lens equation (11) so that the unit distance is te2 



u = z — j-- (17) 

Z\ - a 2 z 2 z 2 



where 



r 



2 

E21 



a 2 = -^p = 1 - d 



' E2 



Case 2): Renormalize the lens equation so that the unit length is Tex- 

1 e- 1 

u = z-- — r — — — (19) 

Z\ — d\Z 2 z 2 



where 



~ ai = r JfL = LA (20) 

r El 6 
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3. Large Separation DSTP Lenses: I >> 1 
3.1. When the Perturber is the First Scatterer 

3.1.1. The Lens Equation 

Let the separation be denoted by i = \x\ — a^i- The coordinate system can be chosen 
such that lens 2 is at the origin, 22 = 0, and lens 1 is on the positive side of the real axis, 



x\ = I. Since it is assumed that £ » 1, the lens equation (17) can be expanded in power 



series in I 1 , assuming that e is not bigger than 0(1), to obtain the following. 

e 1 e _ e 5 2 , Q1 v 

"-£ = 2 -i + ^-^7 (21) 

The lens is made of a point mass (oc 1/z), a constant shear (ex z), and a mass-antimass 
distribution (oc 1/z); the source is shifted by e/l as is the case with the wide separation 
binary lens. (See appendix.) Consider the RHS minus the LHS as a vector field. It is a 
vector field with zeros and poles on the two sphere, and the index of the vector field at 
z ~ 00 results in n + — n_ =0 where n + and n_ are the number of positive and negative 
images respectively. (See RB10.) Thus, the number of images is even and the number of 
negative images is the same as the number of positive images. There are two images for 
uj = 00, namely z = and 00, hence there are two or four images where the latter occurs 
inside the finite size caustic. The finite size caustic occurs because the degeneracy of the 
point caustic of the single lens is broken by the perturbation of the mass M\. The size of 
the caustic curve is calculated below using second order approximation in 



3.1.2. The Critical Curve and Caustic Curve 
The Jacobian of the lens equation is given as 

Jacobian : (22) 

U // 

f = du = l + ff 2 ; g = 8u = ^ + ± (23) 
where d = d/dz and d = d/dz. The Jacobian determinant is 

</ = m 2 -M 2 (24) 

and the eigenvalues of the Jacobian are 
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\± = f R ±w-m 1/2 (25) 

where fu and // are the real and imaginary parts of /. On the critical curve, where J = 0, 
one or both of the eigenvalues are zero because J is the product of the eigenvalues. 

A± = f R ± f R (26) 

Thus A_ vanishes on the critical curve, and A + also vanishes if f R = 0. Note that / = f R = 1 
in the case of the binary lens, and A + never vanishes. Here fn>0 because i » 1. 

The lens system is simple enough so that the critical curve can be explicitly written out 
as a simple function. If we set z = re te , the critical condition is given by the following in the 
linear approximation in e/£ 2 . 

r = l + |^cos20 (27) 

Compared to the circular critical curve r = 1 of the main (single) lens, the critical curve is 
slightly squeezed in a quadrupolar manner. Note that every point of the entire ring r = 1 



of the single lens is a precusp (i.e., mapped to a cusp). The curve in eq.(27) is circular 

[drjdQ = 0) at four points: 0, 7r/2, tt, and 3n/2, and they are expected to be the precusps. 

It is the indeed the case as will be shown shortly. The size of the caustic can be estimated 

by calculating the cusp positions using the lens equation. The precusps along the real axis, 

6 = and tt, are mapped to cusp points on the real axis, and the length of the caustic along 

the real axis is obtained as the difference between the cusp positions. 

4ed ed 
Auj real = u(x) - lo(-x) = — ; x = 1 + — (28) 

The length of the caustic in the direction parallel to the imaginary axis is given as the 
absolute value of the following. 

4ed ed 
4 Mi 2 

|Au; rea /| = \Auji ma g\ = £2~^~" (30) 

Thus the quadroid caustic is equilateral and the orientation of the caustic is opposite to the 
critical curve. In comparison to the wide separation binary lens, for which d = 1, the size 
of the caustic is smaller by factor d 2 . See the Appendix for the wide separation binary. If 
the lens elements are evenly distributed in distance between the observer and the source, 
then d = 1/4, and the caustic shrinks by 1/16. It is substantial, and it demonstrates that 
perturbation of microlensing events by a "rogue" mass in an angular proximity should be 
estimated by using a proper double scattering lens equation. It has been the practice that 
all possible perturbers are universally thrown into constant shear corrections, or constant 
shear and constant convergence. 
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3.1.3. Cusps 

On the critical curve, > 0, hence A_ = is responsible for J = 0. Thus, if e + and e_ 
are the eigenvectors coresponding to A + and A_ respectivley, then e_ is the critical direction. 
If we consider drawing the caustic curve by mapping the critical curve by the lens equation, 
only the non-critical (e+) component of the tangent vector of the critical curve is mapped to 
the tangent of the caustic curve because of the criticality condition. Thus the caustic curve 



is tangent to the eigendirection of A + (Rhie 1999, 2001). If the tangent to the critical curve 



is parallel to the critical direction, the tangent mapped to the caustic curve is zero and the 
progression of the caustic curve stops and forms a cusp. In the next moment, the non-critical 
component is picked up and the caustic curve turns around changing the direction by ir. If 
p is the parameter of the critical curve, the tangent to the curve is determined by 

Q= dJ = dz ±dJ+ dz. d _ J (31) 
dp dp dp 

where dz± are the increments in the ± eigendirections. The cusp forms when dz + /dp = 0, 
hence = <9_ J. The (unnormalized) eigenvectors are 

e+ocfM; e-ocfM, (32) 



where we can choose u± and v± as 



u + = g; v+ = -if I + (\g\i-fif* 

u-=g; V- = -ifr-(\g\Z-ff) 1/2 (33) 

The Jacobian matrix can be diagnolized using A constructed from the eigenvectors compo- 
nents, 

A= ( U v + l~Y- A' 1 (Jacobian) A = (q + a °) 5 (34) 
and the eigendirection differentials can be written as 

dz+ \ . i f dz A 

constant A I (35) 



dz_ J \ dz r 

where constant is a real constant. Thence d-Z = constant w_ and and d-Z = constant i>_, 
and the cusps are found from the cusp condition. 

= 9_J = d-zdJ + d^zdJ = u^dJ + v_dJ (36) 

Straightforward calculations show that, in the second order in i' 1 , d-J = for 8 = 0, tt/2, 



it, and 37r/2 of the critical curve in eq.(27). Therefore they are precusps in the second order 
approximation. 
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3.1.4- Images 



Set z = re and £ = lo — e/£. The lens equation for the shifted source position is given 



by 



( = + r]e 



where £ and 77 are functions of r. 



If we let C = Ci + 



& 1 e ( 3i 
£ = r --; 77= — r 



Ci = (f + ^) cos 0; C2 = (C - »7) sin 9, 



and an equation for r is obained. 



(37) 



(38) 



(39) 



(40) 



There are two or four solutions to the equation (40), which indicates that there are two 



solutions outside the quadroid caustic and four solutions inside the caustic. When the source 
is inside the caustic, the images are all at r ~ 1. They are the four bright images that form 
around the critical curve r ~ 1, two outside the critical curve in the area of the "squeezed" 
and two inisde the critical curve in the area of the "bulged". For example, ( = is inside 
the caustic, and the four images are on the real axis and the imaginary axis. 

* = fc = - 
2 2 



1,2 



i + ^(i-a a ) 



(41) 



3,4 



a 2 



9, = 0, 9 4 



71 



(42) 



The radius is bigger than r 3i4 because a 2 < 1 for a double scattering lens with d < 1. 
Generally, the radii of the images are different. Figure [3] shows a case: £ = 100, d = 1/3, 
and Mi = M 2 ; d = £ 2 = 2. x 10~ 5 . The angle for the radius r of each image is determined 



from eq.(39). The four images of a finite size source filling the caustic form more or less a 
circular ring with finite thickness threaded by r = 1. 



3.2. When the Perturber is the Last Scatterer 



3.2.1. The Lens Equation 



Power-expand eq.(19) in i assuming \z\ ~ 1 because we are interested in the region 
around the critical curve. Keep the terms up to the second order in £~ l because the size 
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of the caustic is of the second order. It will be shown that the linear order perturbation 
shifts the position of the caustic, by d(e£)~ 1 , but does not break the degeneracy of the point 
caustic. 



1 1 5i a\Z 



7,2 



£ 2 z 3 + e£ 2 



(43) 



The lens is made of a point mass, a constant shear, and a whole variety of multipoles. From 
the index of the vector field with zeros and poles, it is obtained that n + — ri- = —2 where 
n± is the number of positive/negative images. Thus the number of images is even. For 
ijj = oo, there are four images, one at z = oo and three degenerate images at z = 0. As w 
moves toward the lenses, the three degenerate images individualizes. It is expected that the 
number of images is four outside the caustic and six inside so that n + — n_ = —2. It is known 
that the original lens equation (eq.(17) or (19)) without approximations produces four or 
six images (Erdl and Schneider 1993 Petters 1997), and RB10 succeeded in deriving the 
sixth order analytic polynomial equation from the lens equaiton. However, the approximate 
lens equation (43) has a third order pole, and we will see that two images are dark images 
remaining near the pole. We refer to them as ignorable images. 



3.2.2. The Critical Curve and Caustic 

Jacobian matrix components are 

„ , a x = 1 2oi 2a 1 z 3a? 1 

f = duJ = 1 + ¥¥ ; 9 = duJ = ¥-Jz^-¥¥ + ¥¥ + d- 2 - (44) 

Set z = re %e and the Jacobian determinant J = \f\ 2 — \g\ 2 can be computed up to the second 
order. 

1 4ai cos 9 2aicos6> /l 2\ 2cos6' Aa\ Qa\cos29 

~ ^4 ^ ' ~jfl I ^2 ^4 J £ £2 r 2 £2 r 6 £2~^6 v > 

In the linear order, 

If 4(2! cos #\ 
Jlinear = 1-^1 ^ J , (46) 

and the critical curve, Ju nea r = 0, is given by 

4a i cos# c . n . ._. 

— — = r - r 5 ; r^O. (47) 

Using graph of the RHS and the fact that | cos^l < 1, it is found that the solution space is 
near r = and r — 1. We are interested in the critical curve with r ~ 1. Let r = 1 — 5 and 
obtain the critical curve in the linear order in 5. 

ai cos 9 



(48 



13 



It is a cardioid even though it is hard to distinguish from a circle because of the small 
coefficient of the cos 9 term. The whole critical curve is mapped by the approximate lens 



equation (43) to a point u> = d(e£)~ l , hence the caustic is a point caustic in the linear order. 
It is shifted from that of the single lens. The position is different by factor d from the center 
of the caustic of the large separation binary lens which is (e£) _1 . 

We need the second order, and the second order depends on cos 29. The combination 
of cos 9 terms and cos 2(9 terms produces an "egg-shape" curve that resembles the familiar 
quaroid critical curve of the binary lens. Set r = 1 — S and 

*=J + | (*») 



and find A and B from the full Jacobian determinant in eq.(45). 



A = 5icos0; B = ^cos2#- ^cos2fl+ ^(1 -cos2fl) (50) 

The critical curve r(9) is a linear function of cos# and cos 29, hence it has the shape of 
an asymmetric peanut (squeezed in at 9 = and it) or a pear depending on whether the 
coefficient of cos 2$ is negative or positive. (The coefficient of cos# is negative.) 

For 9 = and ir, dr/d9 = 0, and they are suspected to be precusps. They are, as can 
be shown by calculating <9_J as was done for case 1). The other two precusps can be seen 
numerically to occur not exactly but practically at 9 = tt /2 and 37r/2. The cusps on the real 
axis will be used to estimate the size of the caustic. 

a Ad AM 2 

Au Teal = ^-^ = - = -— (51) 

It is of order l/£ 2 and is the same size as the large separation binary lens. But the shape of 
the critical curve is different from that of a binary lens as was mentioned above. Such a nice 
result should have a physical interpretation which escapes our mind currenlty. It should be 
worth pondering in an idle time. 

The majority of the microlensing toward the Galactic bulge is bulge-bulge lensing and 
the microlensing event can be perturbed by a foregroud star. The highest magnification 
microlensing event observed to date is of the total magnification 2400, which means that 
the impact distance is 4.17 x 10~ 4 Einstein ring radius. If M\ = M 2 and I = 100, then the 
"radius" of the caustic is 2 x 10~ 4 . If the main lens has solar mass and is at Dn = 3D s /4, 
the Einstein ring radius is tei = 2310sec. If the source star is a sun-like star with the solar 
radius (2.32 sec), it is 10~ 3 in units of the Einstein ring radius. The caustic is completely 
inside the source star radius. Thus perturbers with £ < 1/100 are expected to contribute to 
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measurable effects. Of course, a more massive perturber (large M 2 ) perturbs more strongly, 
and when it is much larger than the main lens, the calculations done here are not proper 
because we would need higher order terms. 



3.2.3. Cusps 



In order to see that the critical points given by eq.(45) are precusps where 9 
calculate cL J for 9 



and 7r, 



and 7r. d- J = 
is easy to see that for 9 = and n, v 



u_dJ + V-dJ where u_ and i>_ are given in eq.(33). It 
_ = — \g\ = —g = iz_, and d J = 8 J in the second order 



in l/£. Therefore, cL J = 0, and the points with 9 = and it are precusps. 



Images and Ignorables 
In the linear order, the lens equation reads as follows. 



1 1 &i 



It produces odd number of images with one more negative image than the positive image as 
one can see from the pole at z = oo and a double pole at z = 0: n + — n_ = — 1. For c<j = cxo, 
there are three images, one at z = oo and two degenerate images at z = 0. Let's look at the 
images at u = 1/ (e£) which can be solved easily. As was discussed before, the point caustic 
is at lj = d(e£)~ 1 , and to = l/(e£) is outside the point caustic and generates three images. 

Z = T' X -2£' - X -2£ (53) 

The second image is the positive image located outside the critical curve where it is squeezed 
in. The third image is a negative image located inside the critical curve where it is bulged 
out. The first image is the dim image located near the lens position z = 0. The first and 
third images become degenerate when u — > oo. One can see that the third image moves 
fast and the first image slowly. Because the caustic is a point caustic, all the three image 
trajectories are continuous. Physically, the first image is ignorable, and practically there are 
two images as is the case in the single lens. 



The full lens equation (in the second order) in eq.(43) has a triple pole at the main lens 
position z = which are image positions of uj = oo. Two of them are ignorable because 
they are confined to the very close proximity to the lens position and have negligible fluxes 
for finite uj. In order to find the approximate positions of the ignorable images for small to 



- 15 - 



(in the near proximity of the caustic), set z = Ajl in the lens equation (43) and find A for 
which the large terms (depending on the positive power of £) add to zero. There are two 
solutions. 

z=°±(l±iV3) (54) 

The position of the corresponding source is O {!/£). Therefore there are two or four images 
in practice. 



The lens equation (43) can be converted into an analytic equation and it is a fifth order 
equation when truncated in the second order in l/£. One of the solutions is an ignorable 
image. The caustic is centered at uj ~ d(e£)^ 1 . 



4. The Double Scattering Two Distributed Mass Lens and Large Separation 

Approximation 

A galaxy lens (of finite mass and finite extension) can be expressed in terms of its 
projected mass density Ma where M is the total mass and a is the normalized projected 
mass density. 

a = 1 (55) 



In the case of a point mass, a is the (Dirac) delta function. If M\0\ and M 2 o"2 are the mass 
densities of the galaxy lenses 1 and 2, the double scattering two distributed mass (DSTD) 
lens equation is written as follows where the unit distance is given by te the Einstein ring 
radius of the total (effective) mass. 

f dPx'aAx') f d 2 x'a2{x') 

u = z ~ ei I fe'V 2 (/) ~ 62 J (56) 



Z\ — x' — a 



z 2 — x" 



where Zj = Z Xj '.J — 1,2 clS before and Xj is the center of mass position of the j-th galaxy 
lens. The density functions Oj are real valued functions and d 2 x' is the real 2-d volume 
element: d 2 x' = dx' A dx'/(—2i). 

If galaxy 1 is the perturber, we can set X\ = £ and x 2 = 0, and assume that £ » 1. 
Again here we are interested in the neighborhood of the critical curve where the bright (or 
detectable) images of a quasar or a galaxy of a high redshift are found. Since the perturbed 
critical curve would be a small modification of the critical curve of the unperturbed lens, 
we renormalize the lens equation so that the critical curve would be given by \z\ ~ 1. 
Let's for conenience denote the second deflection angle (multiplied by a distance factor and 
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normalized) by n 2 . 

„ / \ f d 2 x'a 2 (x') . 
n 2 {z-x 2 = / - _ 21 ^ 57 



The equation (56) can be rewritten as follows where the unit distance is given by rg 2 . 

/d 2 x' 'ai(x') , , . 

z — £ — x' — a 2 n 2 (-2; 0) 

where x 2 = and the double scattering parameter a 2 is defined exactly the same way as in 
the case of the point mass lenses but with the distances of the center of the masses. Now 
the critical curve is in the neighborhood of \z\ = 1, and Il 2 can be assumed to be 0(1). 



The equation (58) can be power-expanded in l/£ to obtain 



"-2 = *- o) + - ^n 2 (z ; o) (59) 

Since the dipole moment of the mass distribution a% is zero, the perturbing galaxy 1 behaves 
as a point mass lens. The last term of eq.(59) depends on the double scattering parameter 
a 2 and vanishes for a perturbing galaxy at the same distance as the main lens leaving only 
the constant shear term. It would be an error to ignore the last term if the (first scattering) 
perturbing galaxy is at a different distance. 

If the perturbing galaxy is lens 2, we can set x\ = and x 2 = i. IT 2 (,2;^) is small and 
can be power-expanded in l/£. 

n 2 (z;£) = -( 1 - + ^)+0(l/P) (60) 



Because of the vanishing dipole moment, the perturbing lens behaves as a point mass. The 
DSTD lens equation (56) can be renormalized so that the unit distance is given by tei and 
the critical curve of the unperturbed lens 1 is given by \z\ = 1. 

u = z- V ; - -n 2 (z; £) 61 

J z - x' - a x Tl 2 {z-£) e 

where x 2 = £. This can be power-expanded with Il 2 as the small quantity assuming that the 
mass distribution of galaxy 1 is well confined inside the critical curve. 

- - a = * ~ n " + + (j + -p) n ' 2 " ¥ Ua (62) 

where 

f (Px'a^x') 

Uln = J Jz^W (63) 
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IIii is the deflection due to the galaxy mass 1 and the third term is the constant shear due 
to the perturber. There are also two other terms, proportional to Il 12 and n 13 respectively, 
that depend on the double scattering parameter a\. It would be an error to ignore these last 
two terms if the (last scattering) perturbing galaxy is at a different distance. 

The caustic of the main galaxy is usually of a finite size unless its projected mass 
distribution is circularly symmetric. The perturbation by another galaxy will change the 
shape and size of the caustic, and it is difficult to handle it algebraically in general. We leave 
the perturbations of finite size caustics for future work. 



A. The Binary Lens with £ » 1 

With d = 1, e = \ 1 1 / \ I > . and the binary lens equation is obtained in which the main 
lens is M 2 . 

a^-^-U^f + l.-I+Or') (Al) 

where the lens positions are x 2 = and x± = £ > 0. The resulting lens is made of a 
point mass (ex z~ l ) and a constant shear (oc z); the source is shifted by e/£. The Jacobian 
components are 

/ = 1; 9 = j 2 +j 2 , (A2) 



and its determinant J is 



J = i-^-£ + °<0 (A3) 



where z = re . The critical condition is given by 

2ecos2# 



r 2 -~ 2 - (A4) 
From the simple graph of the RHS, it can be seen that there is a solution space in the 



neighborhood of r = 1. Set r = 1 + 5 and find 5 using eq.(A3) to obtain 



ecos2t9 

^ = 1 + ^2-- ( A5 ) 

It is a peanut-shape curve squeezed along the imaginary axis even though the small coefficient 
e/2£ 2 makes it difficult to discern from a circle. 8 = 0, tt/2, tt, and 3tt/2 are the precusps. 
It can be confirmed by calculating cL J as it was done for an arbitrary d in the main text. 
The cusps are two on the real axis and two on the imaginaray axis. In order to estimate the 
size of the caustic, measure the cusp-to-cusp distances on the real axis and on the imaginary 
axis. 

4e 

Au rea i = U - U n = — (A6) 
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4e 

The quaroid is equilateral and its orientation is opposite to the critical curve. The diagonal 
length of the quadroid will be comapred to that of the large separation DSTP lens caustics. 

4e 4 Mi 

|AUW| - \Au} ima g\ - ~ £2j^ ( A8 ) 
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Source 



Observer 



Fig. 1. — Double scattering lensing is a time-sequential process and there are two cases: 1) 
lens 1 is the perturber; 2) lens 2 is the perturber. 
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Fig. 2. — A diagram of a photon path scattered by lenses 1 and 2 at distances D\\ and 
D12 in sequence. The scattering planes are noncoplanar in general, and the angles should be 
considered as two-vectors in three space. 
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Fig. 3. — The radii of the four images of a source at d = ( 2 = 2 x 10~ 5 . They are all near 1 
- the Einstein ring radius of the main lens M 2 . The lens parameters are £ = 100, d — 1/3, 
and Mi = M 2 . 
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Fig. 4.— M.2 is the main lens at the origin: The "caustic curve" in blue is obtained 
by mapping the approximate critical curve using the approximate lens equation with the 
perturber as the first scatterer. The perturber is to the right on the real axis. The black 
"caustic curve" is from the same approximate critical curve mapped by the exact DSTP lens 
equation. The exact caustic curve of the DSTP lens (not shown) is hardly distinguishable 
from the black curve. The size of the caustic is ~ 0.013 in agreement with the analytic 
formula. 
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Fig. 5.— Mi is the main lens at the origin: The "caustic curve" in blue is obtained 
by mapping the approximate critical curve using the approximate lens equation with the 
perturber as the last scatterer. The perturber is to the right on the real axis. The black 
"caustic curve" is from the same approximate critical curve mapped by the exact DSTP lens 
equation. The exact caustic caustic curve of the DSTP lens is shown in thick red. The slight 
difference of the black curve from the red curve can be discerned in the area around the real 
axis toward the perturber. The "caustic" sizes are ~ 0.02 in agreement with the analytic 
formula. 
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Fig. 6. — The "caustic curve" in blue is the same as in figj5j The exact caustic curve of the 
approximate lens equation for the last scattering perturber is shown in red. 



